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0 Method to improve directional suivey accuracy. 



0 A method is proposed for better determining the azimuth and inclination of a borehole in which a priori 
information regarding the magnitude of the earth's gravitational field strength, its magnetic field strength and 
direction are utilized to constrain the measured components of the earth's gravitational and magnetic field 
vectors. The constrained fit method adjusts the accelerometer and the magnetometer tri-axis measurements to 
optimize agreement between the measured data and the a priori magnetic and gravity field strengths and the 
magnetic dip angle. The adjustment that is selected is that which produces the least increase in the statistical 
measure x^. 
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METHOD TO IMPROVE DIRECTIONAL SURVEY ACCURACY 



During the process of drilling a borehole in the search for producible hydrocarbons, the practice of 
drilling at an angle to the vertical (inclination) and in a specific direction (azimuth) is becoming more and 
more common This so called directional drilling is usually performed in order to target a specific region in 
the earth (ie a specific geologjcal formation) or to avoid hitting an object such as another well. In order to 

5 know where the borehole is being drilled so that its course may he corrected, downhole surveying 
instruments have been developed and employed. Such surveying instruments typically include a multi-axis 
accelerometer and a multi-axis magnetometer which measure at least two and preferably three components 
of the earth's gravitational and magnetic fields respectively. The azimuth and the inclination of the borehole 
at a particular point is then derived in a known manner from the information measured by the accelerometer 

70 and the magnetometer. In this process, it is assumekj that the earth's gravitational and magnetic fields are 
unaffected by stray or spurious fields. If. however, the fields at the location of the measurement have been 
perturbed by extraneous fields, the determination of the azimuth and inclination of the borehole will be 
incorrect so that the driller will not know with the requisite degree of precision where the borehole actually 
is relative to his "target". 

75 US Patent 4,163,324 describes a technique for correcting for errors that may affect a magnetic 
surveying instrument which resuft from stray magnetic fields arising from the equipment in the borehole. In 
that patent it is assumed that all interference is caused by magnetic material in the drillstring and is, 
therefore, axial. No means are provided for verifying the validity of this assumption. If the assumption is 
wrong, then the correction made to the azimuth measurement will also be- wrong. Thus, the described 

20 technique would give incorrect results if the magnetic interference did not lie along the longitudinal axis of 
the tool as would be the case where the interference arose from an adjacent magnetic anomaly or from 
magnetized components in the tool having transverse magnetic fields. 

US Patent 4.510.696 as well as the SPEIAOC publication 13476 "Reduction of Nonmagnetic Drill Collar 
Length Through Magnetic Azimuth Correction Technique" by A.W. Russel, point out that accelerometer and 

25 magnetometer measurements are subject to errors. This is especially true for the magnetometer whose 
measurement will be influenced by the presence of stray magnetic fields originating from magnetic 
materials in the vicinity of the magnetometer. The example cited by that patent is the magnetic field 
originating from the presence of the drill pipe which, unless made of non-magnetically permeable material, 
is likely to have a residual magnetic field. A further example is the stray magnetic fields that are present 

30 due to the proximity of other measurement while drilling equipment containing magnetically permeable 
material. 

While it is possible to place the surveying instrument within a non-magnetic drill collar with the 
instrument spaced sufficiently far from the origin of the stray field for the field to be small enough to have 
little or no Influence on the measurements, such non*magnetic drill collars are costly, indeed, it is the object 

05 of the above referenced patent and SP&IADC paper to calculate the minimum length of non-magnetic drill 
collars (and thereby minimize the cost) that can be employed without incurring this interference. 

Unfortunately, stray magnetic fields frequently result from components in the drilling assembly other 
than the drill collars. For example, various components of a downhole drilling motor or other downhole 
equipment may be magnetically permeable. Spacing these items sufficiently far from the surveying 

40 equipment may unacceptably constrain the design of the bottom hole assembly. 

In an approach different from that of performing a calculation for determining the minimum size of non- 
magnetic drill collars. US Patent 4.682.421 describes a technique for detecting and removing constant 
magnetic biases from the transverse and axial magnetometer measurements. The technique requires the 
operator to make multiple surveys at various roll angles at a constant depth in a wed in order to determine 

45 the cross axial component of the magnetic field bias due to magnetization of the drill collar. 

In this technique, the X-axis magnetometer measurement is plotted versus the Y-axis magnetometer 
measurement (the Z-axis measurement lying along the longitudinal axis of the drill string) for each of the roll 
angles at which measurements are made. The resultant is a circle whose displacement from the origin is 
indicative of the horizontal X and Y magnetic biases. The longitudinal axis is then corrected by first 

50 subtracting the biases from the transverse axes and then computing the total magnetic field and the 
measured magnetic dip angle. The magnitude of the vector difference l>etween the measured magnetic field 
and the tabulated magnetic field (obtained from a priori independent data) is then calculated- The vector 
difference between the measured magnetic field and the tabulated magnetic fiekj is then used to obtain a 
corrected longitudinal magnetic field. Once these corrections are obtained they are applied to any individual 
survey to correct for the magnetic bias introduced into the measurement by the magnetic drill collar. 
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Several assumptions are implicit in the technique of US Patent 4.682,421. One assumption is that the 
transverse magnetic bias errors are fixed relative to the tool and are not due to. for instance, the position of 
some components of a mud motor which have additional internal degrees of freedom. A second assumption 
is that the magnetic bias errors are stable over time so that biases determined at the time of the calibration 
s procedure can be used to correct data collected at other times. This assumption may not be entirely 
accurate in view of the shock and vibrations to which the drill string is subjected during the drilling process. 
A third assumption is that transverse bias errors dominate other causes of magnetic mis-measurement. For 
example, magnetometer alignment and scale factor errors are assumed to be small relative to the bias 
errors from stray fields when this may not be the case at all. 
fo If these other error sources are the dominant cause of the resultant measurement error rather than the 
bias errors, the measured data will be degraded since erroneous biases will be subtracted from each sensor 
measurement. Likewise, if the magnetic biases are unstable, the data will also be corrupted. 

US Patent 4.761.S89 also describes a technique for addressing the problem of the effect of stray 
magnetic interference on a surveying device. This patent, as well as some of the other above mentioned 
IS techniques utilize a priori magnetic field magnitude and dip values to improve their results. None of the 
atx>ve mentioned techniques, however, attempts to take advantage of the additional a priori gravitational 
field strength information nor do they take into consideration the measurement uncertainties of the 
magnetometers and accelerometers. 

The above techniques are therefore directed exclusively at magnetic field interferences and ignore 
20 poteotial errors in the accelerometer measurements. Since none of the above prior techniques takes into 
account measurement uncertainties that are design limitations of the magnetometer and accelerometer 
sensors themselves, they are unable to produce optimum results. Finally, none of the above techniques 
treats the three magnetometer axes and the three accelerometer axes on a equal basis. Rather, in those 
techniques, undue importance is given to the magnetometer z axis. 
25 It would therefore be advantageous to have a technique for determining inclination and azimuth of a 
borehole that would beneficially utilize ail of the a priori information available including gravitational field 
amplitude and measurement uncertainties. Desirably such a technique should perform corrections to the 
magnetometer and accelerometer measurements without requiring measurements from more than one 
orientation and one depth in the well bore in order to be self contained without having to assume that an 
30 interfering field is invariant from the beginning to the end of the survey of the borehole. 

Furthermore it would be advantageous to have a technique which weighted the contributions of the 
measurements in accordance with the relative estimated uncertainties of the measurements but which 
otherwise treated all of the measurement axes on an equal basis. It woukj be of additional advantage to 
have a technique which could be varied in a manner such as to identify a faulty magnetometer or 
35 accelerometer axis. And finally, it would t>e of great advantage to have a technique that would be 
responsive not only to stray fields produced by objects within the borehole but which would have equal 
sensitivity to field anomalies originating outside of the borehole so that information relative to their location 
and distance could be obtained. 

The shortcomings of the prior techniques for determining inclination and azimuth values having 
40 minimum errors from stray magnetic or gravitational fields are overcome by the present invention while at 
the same time providing additional significant t>enefits and objectives. Briefly, a trt-axis accelerometer and a 
tri-axis magnetometer carried by a drill string make measurements of the components of the earth's 
gravitational field and the earth's magnetic fields. These outputs are then corrected according to calibration 
factors and then are modified to be consistent with three a priori geophysical measurements which include 
45 the earth's gravitational field intensity, the earth's magnetic field intensity, and the earth's magnetic dip 
angle. In this method, an ensemble of accelerometer (gi) and magnetometer (hi) outputs at each 
measurement location are generated to be consistent with the a priori constraints. Improvement occurs 
since the three constraints cause the ensemble of the six measurements to have only three degrees of 
freedom. It is the reduction of the number of degrees of freedom which Improves the results. From 
50 calculated outputs. Qi and hi. the inclination and azimuth of the borehole may be calculated using 
conventional formulas for inclination and azimuth. 

Thus, a procedure which imposes a three constraint fit is performed by the method of Lagrange 
multipliers which minimizes the function: 

55 
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with respect to Qi. h,. and X,. where Gi and are the measured and corrected (for bias, scale factor, and 

,0 alignment errors) accelerometer and magnetometer components respective(> Go. H,, and cos, are the a 
priori values of the earth's gravity field intensity, magnetic field intensity, and the cosine of the angle 
between the gravity and the magnetic field directions. The X, are th^ Lagrange multipliers which introduce 
the constraint conditions into the minimization. The o^,( g ) and Oh,( g . h ) are estimates of the uncertainties 
, in the Gi and the H, measurements. In the minimization procedure, nine non^near simultaneous equations 

,5 are produced which are then solved numerically (for example, by standard IMSL routines), (agj and Oh,, are 
functions of orientation since they reflect uncertainties in the alignment and scale factor of the mag- 
netometer and accelerometer axes. For this reason, they are shown in the a»3ve expressions as functions 
of g"and F. since these vectors determine the tool orientation.) 

This solution results in values of gs and the h, (as well as values for the liree Lagrange multipliers x,) 

20 which are improved estimates of the accelerometer and magnetometer o-?puts since they have been 
constrained to be consistent with the a priori data. 

Figure 1 is an illustration of a functional block diagram illustrating in a general way the functions and 
steps performed by a general or specific purpose digital computer in the practce of the invention. 

In figure 1 there is illustrated a geological formation 10 which is being dnied by a conventional drilling 

25 procedure to form a borehole 12. As part of the drill string there is shown a drill collar 11 having therein 
surveying instrumentation which includes a tri-axial accelerometer 14 and a ai-axial magnetometer 16 for 
making measurements of the components of the earth's gravitational and m^netic fields. The outputs of 
the magnetometer and accelerometer are delivered to a downhole processor ;8 which performs calibration 
corrections with respect to bias, scale factor and alignment errors that have previously been determined for 

30 that particular surveying instrument 

Upon undergoing the above corrections, the signals H,. Q,. representng the outputs of the three 
accelerometers and the magnetometers are either further processed down-hoie to obtain determinations of 
inclination and azimuth or are sent up-hole by a mud pulse telemetry system 20 for further processing at 
the surface as illustrated in processor 22. Processor 22 comprises any standard, suitably programmed 

35 special or general purpose digital computer, as for example the POP 1 1.'35 dig.-al computer. 

In the inventive procedure, the previously corrected values of the components of the measured 
magnetic and gravitational fields. G, and H, are then modified at functional bfeck 32 to be consistent with 
three a priori geophysical measurements 28 which include the scalar magnituae of the earth's gravitational 
field (Go) the scalar magnitude of the earth's magnetic field intensity (Ho), and the complement of the 

40 earth's magnetic dip angle (,). These three a priori quantities may be obtained by reference to standard 
look up tables by knowing the latitude and longitude of the location of the well, v by actual measurement at 
the well site. 

The procedure practiced at 32 involves a constrained minimization of the. r function 

Where the and <»^., represent the gravitational and magnetic field uncertainties determined at functional 

'''°'^The^ preferred method of performing the constrained minimization is =y the method of Lagrange 
multipliers which serves to introduce the three a priori constraints (the scaar magnitude of the earth s 
gravitational field (G„). the scalar magnitude of the earth's magnetic field intensty (Ho), and the coniplement 
• of the earth's magnetic dip angle {,) into the minimization. The function is ftus modified as is known in 
the practice of the Lagrange multiplier method to appear as follows: 
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In the above relationships, the ag., and the cn., are the uncertainties of the gravitational and magnetic field 
vectors attributable to the uncertainties of each of the measurennent axes of the acceleronneter and the 
magnetometer. The uncertainties in the measurement axes arise from uncertainties in the bias, scale factor, 
and alignment which are values available from the vendors of the magnetometer and accelerometer 
instrumentation. Additionally, random uncertainty due to the quantization introduced by the digitization of 
the sensor outputs and the sensor*to-drill collar misalignment are included in the sigmas. The gravitational 
and magnetic field vector uncertainties may be derived at functional block 30 from the following relation- 
ships: 



where the yj are sources of uncertainty in the output of the jth accelerometer output (for example, the bias, 
scale factor, or alignment uncertainty), ^is an estimate of the variance of that source of uncertainty, and 
ly\i is an estimate of the covariance between the various sources of uncertainty. Similarly, the u. ol ,. and al 

are analogous expressions for the sources of error affecting the magnetometer outputs and estimates of 
their magnitudes. The measured values of Gj and Hi can be used to evaluate the above expressions for og 
and ah .] with negligible impact oh the minimization of x^. 

As can be seen from the above relationships, the uncertainties of the magnetic and gravitational field 
components are dependent on the magnitude of the measured components Hj and Gj and therefore must 
be redetermined for each of the surveys performed in the txDrehole. As a result of this and other 
computational complexities, the preferred method of practicing the invention is in a computer 22 located at 
the earth's surface. However, downhole processing in processor 18. is not to be precluded and might very 
well be the preferred mode were processor 18 to possess sufficient processor and memory capacity. 

In the minimization procedure performed at block 32. the partial derivative of the modified function is 
then taken with respect to each of the gravitational and magnetic field components gi and hj as well as with 
respect to the Lagrange multipliers. K Each of these partial derivatives are then set equal to zero to. obtain 
a set of nine simultaneous equations. The resultant nine non-llnear simultaneous equations are then solved 
(for example,, by a standard numerical routine such as that known to the industry as "IMSL") to obtain 
values of gi and the hj (as well as values for the three Lagrange multipliers x,) which are improved estimates 
of the accelerometer and magnetometer outputs gj and hi. 

Since, in this procedure, the results g; and hi have been constrained to be consistent with the a priori 
data, the degrees of freedom have been reduced from six to three. As a result only three unknowns may 
be determined by this technique, only two of which can be from either the magnetometer or the 
accelerometer. Once improved components of gj and hj have been determined, they are used at 34 in 
standard equations to calculate the inclination and azimuth of the borehole. The final inclination and azimuth 
results are output by the processor 22 and recorded by recorder 24 in "log" form. 

The effect of a known interference, such as the proximity of a magnetic drill collar, may be incorporated 
into the above described procedure to produce results for which the perturbation is reduced. Specifically, if 
a magnetic drill collar has a non-zero component (at the location of the magnetometer) extending along the 
longitudinal axis of the tool 1l. (ie the Z axis of the magnetometer), the magnitude of the interfering field 
may be approximated and used to increase the tool's bias uncertainty at 26. It has been determined 
through tests and modeling that the above described procedure is not very sensitive to the exact value of 
the interfering field so that a factor of two approximation will generally suffice. Furthermore, it has been 
discovered that, where the surveying tool is housed in a non-magnetic drill collar, the effects of spurious 



EP 0 384 537 A1 



magnetic fields from near-by magnetic drill collars in the directions normal to the longitudinal axis (the X 
- and Y axes) of the tool 1 1 are so slight that they may he ignored with confidence 

M »t! ^ ""^^^ °' the precision of the results, and therefore of the accuracy of the initial input data G, and 
H. the results may be substituted back into the equation above to calculate a x= value. Where the 

the ,n,t.al data .s so inaccurate that It should not be considered to be a reliable survey and possibly 
discarded Such .naccurac.es may arise in a number of ways such as by the movement of the surveying 
instrumentation dunng the process of measurement or the proximity of a variable source of magnetic 
interference such as a rotating component of a mud motor. 

hinh'llill^nf 'i' V ?f H ^""""^ °' survey instrument axes has failed resulting in an 

high value of x^. If a failed axis .s suspected, the process may be modified as follows to identify the faulty 
axis, which IS subsequently not used. The quantity is repeatedly reformulated and the minimization 
process repeated, with an artificially large uncertainty, o. substituted for the actual uncertainty of one of the 
SIX possible axes. Providing an artificially high uncertainty for a specific axis has the effect of minimizing the 
rs contrrbution of the data from that axis to the value of the quantity. If no appreciable decrease in the 
value IS obtained, it may be concluded that that axis is not the faulty axis and the process is repeated by 
^ providing an artificially high uncertainty for another axis. If. however, the value of x= were to drop 
dramatically, .t .s apparent that the improvement is caused by the lesser contribution from that axis and the 
conclusion is that the faulty axis has been identified. 
20 An additional variation that is available using the above technique is to derive additional information on 
spunous magnetic and gravitational field sources as follows. If it is suspected that the borehole is nearing 
an adjacent well in which there is a magnetically permeable casing or drill pipe and that the borehole is not 
changing its onentation. the dip and azimuth values may be used as a portion of the a priori data In this 
manner, the results of a number of surveys may be compared with one another. The variations in the 
2S results may be assumed to be attributable to the changing proximity of the anomaly as the borehole 
conta.n.ng the surveying equipment changes its position relative to the anomaly. Repeated application of 
this technique may enable determination of the direction and possibly distance to a nearby magnetic field 
source. For example, after applying this technique, the differences between the measured and fit values of 
the magnetic field vector at several known distances along the well bore can be used to determine the 
30 direction and pole strength of the source of the magnetic anomaly. 

While preferred embodiments have been shown and described, viarious modifications and substitutions 
may be made thereto without departing from the spirit and scope of the invention. Accordingly it is to be 
understood that the present invention has been described by way of illustrations and not limitation 
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1. A method for surveying a borehole formed through subsurface geological earth formations compris- 
ing the steps of: 

a. at an in situ location in said borehole, measuring a plurality of components of the gravitational field 
strength and a plurality of components of the magnetic field strength: 

b. determining an a priori value for the earth's gravitational field at the latitude and longitude of the 
borehole: and 

c. in response to said measured components and to said a priori value, determining the azimuth and 
«. the .nclination of the borehole at said in situ location in the borehole. 

2. A method for surveying a borehole formed through subsurface geological earth formations compris- 
ing the steps of: 

a. at an in situ location in said borehole, measuring, with magnetic and gravitational field responsive 
measuring instruments, a plurality of components of the gravitational field strength and a plurality of 

so components of the magnetic field strength; 

b. determining measurement uncertainties for each of said gravitational and magnetic field comoo- 
nents: and 

c. in response to said measured components and to said measurement uncertainties, determining the 
azimuth and the inclination of the borehole at said in situ location in the borehole. 

3. The method as recited in claim 2 further including the step of determining an a priori value for the 
earth's gravitational field strength and wherein said step of determining the azimuth and the inclination of 
the borehole at said in situ location is responsive to said a priori gravitational fiekJ strength. 

4. The method as recited in one of claims 1 or 3 further including the step of determining a priori values 

6 



55 



EP 0 384 537 A1 



for the earth's magnetic field strength and wherein said step of determining the azimuth of the borehole at 
said in situ location is responsive to said a priori magnetic field values. 

5. The method as recited in claim 1 further including the step of determining the measurement 
uncertainties for each of said gravitational and magnetic field components and wherein said step of 

5 determining the azimuth and the inclination of the borehole at said in situ location is responsive to said 
measurement uncertainties. 

6. The method as recited in one of claims 2 or 5 wherein the measurement uncertainty of each of said 
magnetic and gravitational field components is determined from the measurement uncertainties of said 
magnetic and gravitational field responsive measuring instruments including scale factor, bias, and align- 

ro ment of the measuring sensors from which said gravitational and magnetic components are obtained. 

7. The method as recited in claim 4 further including determining a value representative of the dot 
product between the earth's gravitational and magnetic field and wherein said step of determining the 
azimuth and the Inclination of the borehole at said in situ location Is responsive to said value representative 
of said dot product. 

fS 8. The method as recited in claim 4 wherein said step of determining the azimuth and the inclination of 
the borehole at said in situ location includes the step of imposing a constrained minimization subject to 
constraints imposed by one of the memt>ers of the group comprising said a priori gravitational field value, 
said magnetic field value and their dot product. 

9. The method as recited in claim 8 in which said constrained minimization treats each of said 
20 constraints and said magnetic and gravitational field components symmetncaliy. 

10. The method as recited in claim 6 wherein each of said magnetic and gravitational field components 
are weighted in accordance with its respective measurement uncertainty. 

1 1 . The method as recited in claim 4 wherein the quantity 
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is minimized subject to the constraints 

iml 

35 S 



45 



so 



where 

Gi = the measured components of the gravitational field 

Hi = the measured components of the magnetic field 

gj - the improved values of the components of the gravitational Held 

hi ~ the improved values of the components of the magnetic field 

Go the a priori magnitude of the gravitational field 

Mo = the a priori magnitude of the magnetic field 

ii = the a priori inclination of the magnetic field 

<7q.i = the uncertainties of the components of the gravitational field 

9h.t = the uncertainties of the components of the magnetic field. 

12. The method as recited in claim 11 wherein said- minimization includes the step of minimizing the 
distribution 
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with respect to Qu h, and Xj, 

where Xi . X2, and X3 are the Lagrangian multipliers. 

13. The method as recited in claim 11 wherein said minimization is repeated in order to identify an axis 
of said magnetic and gravitational field measuring instruments which is faulty, each minimization repetition 
being performed with an uncertainty on successively different axes which is large relative to the uncertain- 
ties utilized for the other of said axes, whereby the faulty axis is identified when becomes small. 
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